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It is shown that attempts to accurately deduce the magnetic penetration depth A of overdoped 
BaFei.g2Coo.igAs2 single crystals by transverse-field muon spin rotation (TF-/iSR) are thwarted 
by field-induced magnetic order and strong vortex-lattice disorder. We explain how substantial 
deviations from the magnetic field distribution of a nearly perfect vortex lattice by one or both of 
these factors is also significant for other iron-based superconductors, and this introduces considerable 
uncertainty in the values of A obtained by TF-/iSR. 

PACS numbers: 74.72. Bk, 74.25.Ha, 76.75. +i 



TF-/iSR is routinely used to determine the magnetic 
penetration depth A of type-II superconductors in the 
vortex state for the purpose of obtaining indirect infor- 
mation on the energy gap structure The magnetic 
field distribution n{B) in the sample is determined by 
detecting the decay positrons from implanted positive 
muons that locally probe the internal fields, and A is 
subsequently determined by modeling the contribution 
of the vortex lattice (VL) to n(B). However, even in 
conventional superconductors the VL contribution is not 
known a priori, and one must rely on phenomenological 
models to deduce what is really an "effective" penetration 
depth A. One reason for this is that only cumbersome 
microscopic theories account for the effects of low-energy 
excitations on n(B) 0). Extrapolating low-temperature 
measurements of A to zero field to eliminate the effects of 
intervortex transfer of quasiparticles, as well as nonlocal 
and nonlinear effects, has been demonstrated to be an 
accurate way of determining the "true" magnetic pen- 
etration depth A 0, 0]- Yet an underlying assumption 
is always that the VL is highly ordered and that other 
contributions to n(B) are relatively minor. The purpose 
of this Letter is to point out that this is not the case in 
many of the recently discovered iron-based superconduc- 
tors, making a reliable determination of A by TF-/iSR 
extremely difficult. 



Here we report on representative TF-//SR measure- 
ments of BaFei.82Coo.isAs2 (T c — 21 K) single crys- 
tals grown from a FeAs flux, as described elsewhere [5J. 
High-statistics TF-^iSR spectra of 20 million muon decay 
events were collected in magnetic fields H = 0.02 T to 0.5 
T applied transverse to the initial muon spin polarization 
P(t = Q), and parallel to the c-axis of the crystals. The 
TF-/iSR signal is the time evolution of the muon spin 



polarization, and is related to n(B) as follows 

/•oo 

P(t)= / n(B)exp(i lu Bt)dB, (1) 
Jo 

where 7 M is the muon gyromagnetic ratio. Generally, the 
TF-/jSR signal is fit in the time domain, with the inverse 
Fourier transform or "TF-/iSR line shape" given by 

/•oo 

n(B)= / P(t)exp(-i^Bt)dt, (2) 
Jo 

providing a visual approximation of the internal field dis- 
tribution. The field distribution of a perfectly ordered 
VL is characterized by sharp cutoffs at the minimum and 
maximum values of B{v), and a sharp peak at the saddle- 
point value of B(r) |il|. These features are not observed 
in polycrystalline samples, where the orientation of the 
crystal lattice varies with respect to H, but are observed 
in single crystals when a highly-ordered VL exists and 
other contributions to n(B) are relatively minor. 

We have attempted to fit the TF-//SR spectra to a the- 
oretical polarization function P(t) that has been success- 
fully applied to a wide variety of type-II superconductors, 
and utilized in some of the experiments on iron-based 
superconductors. The spatial variation of the field, from 
which n(B) is calculated, is modeled by the analytical 
Ginzburg-Landau (GL) function [l[ 

B(r) =B (l-b )^ ~— , (3) 

g A G 

where G are the reciprocal lattice vectors of an hexag- 
onal VL, b — B/B c2 is the reduced field, B is the av- 
erage internal magnetic field, Ki(u) is a modified Besscl 
function, and u 2 = 2<fG 2 (l + 6 4 )[1 - 26(1 - bf). As 
explained later, P(t) is multiplied by a Gaussian depo- 
larization function exp(— a 2 t 2 ) to account for the effects 
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FIG. 1: (Color online) TF-/jSR line shape of 
BaFei.82Coo.i 8 As2 at # = 0.5 T and T — 3.9 K (green 
circles), (a) The red curve is the Fourier transform of a fit in 
the time domain assuming B(r) of the ideal VL is given by 
Eq. (|3}. In addition to the indicated values of A and £, the 
fit yields a = 0.265 /js _1 and a PM shift of 8.6 G. (b) Fourier 
transform of a fit in the time domain that assumes the model 
of field-induced AF order described in the main text (red 
curve). Some of the fit parameters are shown in Fig. [3] 

of nuclear dipolar fields and frozen random disorder. We 
stress that the fitting parameters £ and A can deviate 
substantially from the "true" coherence length and mag- 
netic penetration depth if other contributions to n(B) 
are significant. An important feature of Eq. ([3]) is that 
it accounts for the finite size of the vortex cores, by gen- 
erating a "high-field" cutoff in n(B). The GL coherence 
length £ a b ~ 26 A calculated from the upper critical field 
H c2 - 50 T of BaFei. 84 Coo.i 6 As 2 with H || c @, repre- 
sents a lower limit for the vortex core radius fij • The core 
size can be much larger if there are spatially extended 
quasiparticle core states associated with either the exis- 
tence of a second smaller superconducting gap or a 
single anisotropic gap Q . Yet fits of the TF-//SR spectra 
of BaFei.82Coo.i8As2 using Eq. ([3]), show no sensitivity 
to the vortex cores at any field and converge with values 
of £ approaching zero. Fig. [T] shows that even at 0.5 T 
where the vortex density is highest, a high-field cutoff 
is not discernible in the TF-/L/.SR line shape. We next 
discuss reasons for this insensitivity to the vortex cores. 

Magnetism — The effective field B^ experienced by the 
muon is a vector sum of various contributions that may 
be static or fluctuating in time. With correlation times 
generally much longer than the muon life time, the nu- 
clear moments constitute a dense static moment system 
that cause a Gaussian-like depolarization of the TF-/jSR 
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FIG. 2: (Color online) (a) Envelopes of TF-/jSR spectra of 
BaFe1.g2Coo.i8 As2 in the normal state at T = 23K. The 
solid curves are fits to a single exponential relaxation func- 
tion G(t) = exp(-At), yielding A = 0.081 ±0.003 /is -1 and 
A = 0.119±0.003 /its -1 at # = 0.02 T and # = 0.5 T, respec- 
tively, (b) TF-/iSR line shapes of BaFe1.82Coo.i8 As2 in the 
superconducting state at # = 0.02 T. The dashed vertical line 
indicates the magnitude of the applied field #. 



spectrum. Yet as shown in Fig. [2ja), BaFei.82Co .isAs2 
exhibits an exponential depolarization above T c that is 
typical of dilute or fast fluctuating electronic moments 
[9j . The latter interpretation is consistent with the obser- 
vation of a paramagnetic (PM) shift of the average inter- 
nal field (B p ) sensed by the muons below T c . This is ev- 
ident in Fig. [2]Jb), where we show representative Fourier 
transforms of P(t) at # = 0.02 T. Instead of the expected 
diamagnetic shift imposed by the superconducting state, 
(B M ) exceeds H. The magnitude of the PM shift in- 
creases with increasing H and/or decreasing T. 

The occurrence of a PM shift in the superconducting 
state of BaFe2-.rCo. r As2 and SrFe2-a;Co x As2 has been 
reported by others [111 . Il2j , and implies an enhancement 
of (B M ) from magnetic order occupying a large volume of 
the sample. Magnetic order exists in underdoped sam- 
ples at H = [l(| , and is apparently induced in overdoped 
samples by the applied field. Yet the effects of magnetism 
on the line width and functional form of n(B) have not 
been considered. A strong relaxation of the TF-/1SR sig- 
nal occurs even in long-range magnetically ordered sys- 
tems, and with decreasing temperature there must be an 
increased broadening of n(B) associated with the growth 
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FIG. 3: (Color online) Results of fits of the TF-/1SR time spec- 
tra of BaFei.82Coo.i8As2 at H — 0.5 T, assuming the model 
of magnetic order described in the main text. Temperature 
dependence of (a) 1/A 2 , (b) the depolarization rates a and 
A, and (c) the amplitude of the additional AF field in the 
vortex-core region. Also shown in (a) are results of fits with- 
out magnetic order and £ = 5 A(blue circles). 



of the correlation time for spin fluctuations. 

Accounting for such magnetism is non-trivial because 
of the spatially-varying superconducting order parame- 
ter and the likelihood that the field-induced magnetism 
occurs in a nematic phase [l3[ . Even so we have achieved 
excellent fits of the TF-/xSR spectra of BaFei.82Coo.isAs2 
to polarization functions that incorporate enhanced mag- 
netism in the vortex core region (e.g. commensurate 
spin-density wave, ferromagnetism, spin-glass), where su- 
perconductivity is suppressed. Here we describe typical 
results for one form of magnetism: To account for line 
broadening from magnetism at higher temperatures, P(t) 
is multiplied by an exponential depolarization function 
exp(— At), as observed above T c . In addition, enhanced 
magnetic order in the vortex cores is modeled by adding 
the following term to Eq. ^ 



B AF (r) = B AF e' 



iK 



iK-i 



K 



(4) 



The K sum is the reciprocal lattice of an antiferromag- 
netic (AF) square iron sublattice of spacing a = 2.8 A, 
Baf is the field amplitude, £af governs the radial decay 
of the field amplitude from the core center, and r and 



r' are the position vectors for 'up' and 'down' spins, re- 
spectively. This kind of magnetic order has the effect of 
smearing the high-field cutoff, and can even introduce a 



low-field tail in n(B) [14J. As indicated by the value of 



£ in Fig. [TJb), fits to this model are sensitive to the vor- 
tex cores. With decreasing temperature, the magnetism- 
induced relaxation evolves from exponential to Gaussian, 
and the magnetic order in the vortex cores is enhanced. 
Consistent with the results Ref. 12j, fits of TF-/1SR spec- 
tra of overdoped BaFe2- a: Co a ;As2 to a model that does 
not include magnetism and is insensitive to the vortex 
cores (i.e. £ fixed to 5 A) yield an unusual linear temper- 
ature of 1/A 2 immediately below T c , and a saturation of 
A at low T . In contrast, fits assuming magnetic order ex- 
hibit a linear temperature dependence well below T c that 
is suggestive of gap nodes. However, these results simply 
demonstrate the ambiguity in modeling such data. With- 
out knowledge of the precise form of the magnetism, our 
model cannot be deemed rigorously valid. Furthermore, 
as we explain next, VL disorder is a serious concern. 

Disorder — Thus far TF-/xSR has been applied to iron- 
based superconductors under the assumption that one 
is probing a fairly well-ordered hexagonal VL. Yet to 
date this has been observed only in KFe2As2 flij . 
Vortex imaging experiments on the i?FeAs(Oi_^F :c ), 
Ai_j;i? x Fe2As2 and AFe2-^Co s As2 families all show a 
hig hly disordered VL indicative of strong bulk pinning 
16 21]. In Fig. ;4j we show the effect of such dis- 
order on the ideal n(B). Using a radial distribution 
function closely resembling that observed in overdoped 
BaFei.8iCoo.igAs2 [2QJ, we have used molecular dynam- 
ics to simulate n(B) of the disordered VL. Although the 
disordered line shape in Fig. 2{b) is asymmetric, it is 
strongly smeared with a field variation greatly exceeding 
that of the perfect VL. 

Small perturbations of the VL by random pinning can 
be handled by convoluting the ideal theoretical line shape 
with a Gaussian distribution of fields [23]. This causes 
a Gaussian depolarization exp(— a 2 t 2 ) of P(t). But for 
polycrystalline samples, n(B) of the perfect VL is nearly 
symmetric and cannot be isolated from a symmetric dis- 
tribution of disorder. Consequently, VL disorder has not 
been accounted for in TF-/iSR studies of pol ycry stalline 
or powdered iron-based superconductors [23M25I ] . Given 
the severity of disorder in these materials and no knowl- 
edge about how this disorder evolves with temperature 
or doping, the accuracy of information deduced about A 
is questionable. Since disorder of rigid flux lines broaden 
n(B), such studies certainly underestimate A. 

While small perturbations of B(r) by vortex pin- 
ning may be accounted for in measurements on single 
crystals, a Gaussian convolution of the ideal n(B) be- 
comes increa sing ly inadequate as the degree of disorder 
is enhanced [26[. In Fig. HJb) we show that Gaussian 
broadening of the ideal line shape does not precisely re- 
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FIG. 4: (Color online) (a) Radial distribution function (RDF) of BaFei.8iCoo.i9As 2 at H — 0.5 T from Ref. [2(j and of the 
disordered VL shown in the lower right generated by molecular dynamics (MD). Note 5000 vortices were used in the MD 
simulation. The horizontal scale is normalized with respect to the intervortex spacing a = 691 A of the perfect hexagonal 
VL. (b) Theoretical simulations of the TF-/xSR line shape of the perfect VL (black curve) and of the disordered VL (red 
curve) corresponding to the RDF shown in (a). The green curve is the line shape of the perfect VL convoluted by a Gaussian 
distribution of fields, corresponding to a— 1.9 /xs _1 . All three simulations assume A = 2000 Aand £=50 A. (c) Same disordered 
line shape shown in (b) and a Gaussian broadened ideal line shape with A = 1600 A, £ = 100 A, and cr = 1.8 /is -1 . The heights 
of the line shapes in (b) and (c) are normalized with repect to the height n max (_B) of the ideal line shape. 



produce n(B) of the disordered VL. More importantly, 
however, because the large disorder-induced broadening 
smears out the high-field cutoff, the fitting parameters 
A and £ are ambiguous. This is illustrated in Fig. @|c), 
where a nearly identical Gaussian broadened line shape 
is obtained for very different values of these parameters. 
Hence disorder introduces considerable uncertainty even 
in measurements on single crystals [ll|, EH, l27l - [3C] | . 

In summary, the effects of magnetic order and/or ran- 
dom frozen disorder of the VL in iron-based supercon- 
ductors introduce considerable uncertainty in values of A 
obtained by TF-/iSR. Unfortunately, these effects cannot 
be modeled in a reliable way. Compounding the problem 
is a lack of information on how the magnetism and VL 
disorder evolve with temperature. Consequently, caution 
is warranted in drawing conclusions about the anisotropy 
of the superconducting gap in these materials from TF- 
/iSR measurements. 
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